Effect of polydispersity to specific absorption rate 
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The predictions of a newly developed Bloch-Bloembergen alike analytic magnetization model are compared 
to experimental results. The effect of size polydipersity on the specific absorption loss is demonstrated for the 
magnetic nanoparticles containing media. Specific absorption rate shows resonance like behivior as a function 
of particle size. The obtained results are in excellent agreement with experimental data. The dominace of the 
Neel relaxation over the Brownian one is demonstrated. 
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I. INTRODUCTION 

Magnetic fluid hyperthermia (MFH) has gained wide inter- 
est in its applicability in medical sciences Enhancing spe- 
cific absorption loss (SAR) by core-shell nanostructures, see 
e -g- JUS] or via anisotropy J3llt] are also extensively studied. 
Furthermore, efficiency of circularly polarized field opposed 
to linearly polarized field is also investigated by models |@- 
[l2tl and experiments lfl3i [Till as well. An interesting study of 
the chain formation of nanoparticles can be read in the paper 
of He lfl5ll . examining the possibility of data storage. 

In MFH, the basic question is how to describe energy 
losses. In the literature relaxation and hysteresis losses are 
distinguished though both proportional to the area of hystere- 
sis curve. See also the comment about it in the paper of Carrey 
et al. flrjll . in which they also argue the necessity of distinc- 
tion. The relaxation losses could originate either turning the 
single-domain particle with its magnetic momentum, called 
Brownian relaxation, or the nanoparticle itself is fixed but its 
magnetic momentum aligns to the external field, as in the case 
of Neel relaxation. There is a vivid dispute, whether both re- 
laxations contribute to the losses and if so, to what extent. 
Also in the paper of Wang et a/. lfl7ll . they compared SAR val- 
ues for nanoparticles with and without polymerization, and 
they found no change in the SAR value excluding the possi- 
bility of Brownian relaxation. Furthermore from theoretical 
considerations, at typical size of nanoparticles up to few tens 
nanometers of diameters, depending on the experimental con- 
ditions, Neel relaxation is regarded to be the dominant pro- 
cess. 

It is also a question, whether a single independent particle 
or statistical ensemble of particles should be taken into ac- 
count. The single particle picture has the advantage of ab ini- 
tio description, though in this case no temperature and there- 
fore thermal fluctuation is included yet Ja, HI, while in the 
statistical picture the treatment is phenomenological lfl8i[l9ll . 
Shliomis itlal provided a complex equation of motion of the 
magnetization from hydrodynamical considerations, but as we 
saw till now, radical simplifications were used leaving only 
the Debye-term of the equation, see e.g. paper of Cantillon- 
Murphy et al. J3]. Stochasticity can also be included U for a 



more elaborated picture of possible relaxation processes. 

Aim of this paper is to test the validity of magnetization 
dynamics model with Larmor-precession term included ll20ll 
in comparison with experimental results llr7[ I2l"l - l23ll. In the 
model the so-called Bloch-Bloembergen equation luil l25ll is 
rewritten for the single domain nanoparticle magnetization. 
We show that experimental results can be understood by the 
polidispersity of size of nanoparticles. Doe to analytic solu- 
tion for both linearly and circularly polarized fields, the model 
is unique being valid without restrictions to any parameter 
value. 

For example, in the paper of Mehdaoui et aZ. ll2lll . a com- 
bined theoretical and experimental study was conducted, in 
which linear response theory and Stoner-Wohlfarth model 
were used, according to their range of validity. We provided 
reliable SAR values with our model both in the range of linear 
response theory and Stoner-Wohlfarth model. 

For practical reason, we focus on hyperthermia applica- 
tion and carried out our analysis and evaluation at lower field 
strength and frequency region. In this case we can claim the 
equivalence of linearly and circularly polarized field for SAR. 
Finally, dominance of Neel relaxation with respect to Brown- 
ian relaxation is also discussed. 



II. THEORETICAL BACKGROUND 
A. Basic equations of motion 

The basic equation of motion of the magnetization is given 

as, 



dM/dt = 7M x B. 



(1) 
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Here 7 is the gyromagnetic ratio with value 7 = —1.76 x 
10 11 Am 2 /Js. It comes from eq.([TJ that any change of the 
magnetization is perpendicular to M and for constant B 
the angle between the two vectors M and B is also con- 
stant. The solution is the precessing M around B for this 
equation of motion, yielding the Larmor precession, where 
= yM x B/Mj_. Here M± is the projection of M on 
the plane perpendicular to B. The Larmor frequency is de- 
fined as a positive quantity, lul = \"f\B. 

Shliomis in 1974 lfl8ll has suggested the equation of expo- 
nential relaxation to be written for describing the behaviour of 
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magnetic nanoparticles in ferrofluids at given conditions, 



ized fields: 



dM M - M P 



dt 



(2) 



where M is the average magnetization of the particles, 
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V is the particle volume, e# = H/H is the unit vector 
along H and Ms is the saturation magnetization of the colloid 
Ms = 4>Md, with volume fraction and the single-domain 
magnetic nanoparticle magnetization M^, The Langevin 
function, C(x) = coth(x) — 1/x, gives the magnitude of the 
magnetization in thermal equilibrium, where M eq must be an 
ensemble average and so is M. The eq.© is the reduced form 
of the Shliomis relaxation equation 1 1 <S| and often called De- 
bye relaxation equation. Susceptibility is frequently approxi- 
mated with Rosensweig's chord susceptibility l26ll . 
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instead of the Langevin function appearing in eq.©, where 
actual magnetic field is used contrary to its maximum. 

The Debye relaxation equation, eq.©, can be regarded 
as a simplified version of the Bloch-Bloembergen equation 
I24H25I . which describes nuclear magnetic resonance experi- 
ments l24ll and ferromagnetic resonance i25ll . The latter does 
include the gyromagnetic torque, pojM x H. The Bloch- 
Bloembergen alike formalism was used in paper |20|| for the 
description of the dynamics of magnetization of paramagnetic 
nanoparticles in ferrofluids. It is pointed out that at field 
strengths and frequencies used in MFH, the linearly and circu- 
larly polarized fields can be equivalent up to a normalization 
constant, which holds for the interested hyperthermia region 
in this paper. 



B. Model equations 

Analytical solutions of Bloch-Bloembergen alike equations 
are briefly summarized here for circularly and linearly polar- 



dM M-M eq 

— — = fialM x H q -, 

dt t 



(5) 



with M eq as given in eq.((3]l. For details of the derivation of 
equations, see paper i20ll . 

Setting rotating magnetic field, H x = Ho cos(wi); H y = 
H sm(u>t), the solution is: 



M(t)=M eq (H ) 



fi 2 1 + {fir) 2 



'us cos(wt) + fi 2 r sin(u;i) N 
ui sin(wt) — fi 2 r cos(wi) I , 



-UJ L 



(6) 

where fi = ^lo 2 +lo\ and lol = hq^Hq. Calculating the 
energy loss per cycle, 
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The specific absorption rate (SAR) is defined as energy loss 
for unit time and mass: 



SAR 



Eus p M eq (H )H lu 2 fir 
2ir~p ~ p ~fi 1+ (fir) 



(8) 



In case of linearly polarized field, H(t) = 
(0,0,HoCOs(uit)), the solution for M z (t) is found in 
the form of effective magnetization: 



_j cos(wi) + lot sin(wt) 
1 + (ojt) 2 ; 



M?'{t) = 2M S £ 7^2 B ^C 2 

(9) 

where £?2m and ^ 2m_1 are the Bernoulli numbers [Ref. Hjl 
p. 1040 and p. 1045] and the argument of the Langevin function 
C = poHoMdV /kT respectively. 
The energy loss per cycle is 



J 



r to+2ir/u> ^jj 

-p I M ■ —^dt 



to 



dt 



M e Jf ■ -^dt 



KpoXloss(() H 



1 + (wr) 2 



(10) 



r 



and the corresponding SAR is 
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r = 2.8nm 


r = 9.85nm 


r = 13.75nm 


O A 7~1 / T X 7 I \ 

OAR M ehdaoui(W 1 g) 


20 


290 


40 


SAR N {W/g) 


21 


291 


42.5 


u M ehdaoui 


0.2 


0.09 


0.2 


calc 


0.16 


0.13 


0.19 


±8 ca i c (nm) 


-1.2, +1.8 


-3.2, +5.3 


-6.3, 11.8 





Wang et al. 


Hergt et al. 


Suto et al. 


SARpap er {W/g) 


123 


45 


22 


SAR N {W/g) 


3 


19 


30 


SAR tot {W/g) 


3 


19 


35 


SAR B {W/g) 


767 


1313* 


5.5 


r(nm) 


2.5 


2.8 


3.125 



TABLE I: SAR (W/g) values are compared for different nanoparti- 
cles sizes. SAR from paper Mehdaoui et al. l2lll (SARMehdaoui) 
and SAR calculated taking Neel-relaxation time (SARn) only are 
shown. Log-normal distribution is taken. Their standard deviations 
are also given by Mehdaoui et al. (a Mehdaoui) and by us (a catc ). 
The possible deviations from the mean value (±5 ca i c ) are given in 
the last row. 



TABLE II: Without polydispersity, SAR (W/g) values are compared 
for different papers. SAR from Wang et al. u3§\, Hergt et al. I22I1 
and Suto et al. llzll are listed. Samples were for Wang et al. size of 
r = 5nm and for Suto et al. the so-called SampleA. SAR calcu- 
lated taking Neel-relaxation time (SARm) only, combined Neel- and 
Brownian-relaxation time (SARtot) and Brownian-relaxation time 
(SARb) only. * No assumption for the value of hydrostatic volume 
is taken. 



III. COMPARISONS WITH EXPERIMENTAL RESULTS 
AND DISCUSSION 



In this chapter our model predictions are compered to ex- 
perimental data. All the parameters can be found in more 
detail in references. First of all, we highlight the paper of 
Mehdaoui et al. Ir2lll . because they included very detailed 
information on polydispersity of their samples and thorough 
analysis of their measurement. It provided great opportunity 
to test our model with polydispersity built in. Excellent agree- 
ment was found and it encouraged us to analize other resutls 
when polydispersity is mentioned but no additional data are 
given. Log-normal distribution was used like in the paper of 
Mehdaoui et al. I2H1 . as particle size could take only non- 
negative values. Fitting log-normal distribution to experimen- 
tal data, we can conclude not only the possible standard devia- 
tions from the average value but also might distinguish among 
samples fabricated by different processes. 

Analyzing samples from Mehdaoui et al., see Table Q] we 
found excellent agreement both for SAR and standard devi- 
ation reproducing also the particle size distribution with our 
model. Mehdaoui et al. used samples with mean radius 
2.8nm, 9.85nm and 13.75nm and found ferromagnetic behav- 
ior with r > 7nm. Table U shows an other very important 
feature of our model providing reliable SAR values out of 
the range of linear response theory, where originally Stoner- 
Wohlfarth model was used instead. Accordingly we carried 
out analytical simulation for direct comparisons with experi- 
ments instead of numerical ones. Finally, we included poly- 
dispersity in our model. 

In Table HI| experimental SAR data from Wang et al. lfl7ll . 
Hergt et al. l22ll and Suto et al. l23ll are compared to model 
predictions calculated without polydispersity. It can be seen, 
the predictions are at the order of experiments, but inconsis- 
tently under- or overestimating them. Data also show that tak- 
ing into account only Brownian relaxation results in extremely 
poor agreement among model predictions and experiments. 
Including polydispersity in calculations results in drastic im- 
provement in the consistency of experimental data and model 
prediction, the error is less than 5% in most cases. In case of 
pure Brownian relaxation process the agreement is poor again. 
Based on these results, the scenario of pure Brownian relax- 
ation process can be excluded. Whether Brownian relaxation 





Wang et al. 


Hergt et al. 


Suto et al. 


SAR Pap er(W/g) 


123 


45 


22 


SAR N (W/g) 


122 


47 


19 


SAR tot (W/g) 


129 


47 


23 


SAR B (W/g) 


146 


146* 


5.8 


r(nm) 


2.5 


2.8 


3.125 


f calc 


0.27 


0.14 


0.25 


±5 C aic(nm) 


-1.5, +6.5 


-1.0, +3.0 


-2.0, +5.5 



TABLE III: With polydispersity, SAR (W/g) values are compared for 
different papers. SAR from Wang et al. Uj\\, Hergt et al. |22|] and 
Suto et al. l23ll are listed. Samples were for Wang et al. size of 
r = 5nm and for Suto et al. the so-called SampleA. SAR calcu- 
lated taking Neel-relaxation time (SARm) only, combined Neel- and 
Brownian-relaxation time (SAR to t) and Brownian-relaxation time 
(SARb) only. Their standard deviations (a ca i c ) are calculated by 
us only. The possible deviations from the mean value (±S C aic) are 
given in the last row. * No assumption for the value of hydrostatic 
volume is taken. 



process is still contributing to the relaxation processes or only 
Neel relaxation should be considered, cannot be decided. The 
first possibility slightly overestimates, the latter one underes- 
timates the experiments. 

Finally, we compared our model predictions to a series of 
data taken from Wang et al. lfl7ll . We experienced the same 
tendencies as before. The consideration of Brownian relax- 





r = 5nm 


r = 4nm 


r = 3nm 


SAR Wang (W/g) 


123 


78 


50 


SAR N (W/g) 


122 


76 


47 


SAR tot (W/g) 


129 


82 


53 


SAR B (W/g) 


146 


496 


304 


O calc 


0.27 


0.36 


0.47 


±5 C aic(nm) 


-2.6, +7.4 


-2.0, +7.5 


-2.0, +8.5 



TABLE IV: SAR (W/g) values are compared for different nanopar- 
ticles sizes. SAR from paper Wang et al. IT7I1 (SARwang) and 
SAR calculated taking Neel-relaxation time (SARm) only, com- 
bined Neel- and Brownian-relaxation time (SARtot) and Brownian- 
relaxation time (SARb) only are shown. Log-normal distribution is 
taken. Their standard deviations (a ca i c ) are calculated by us only. 
The possible deviations from the mean value (±S ca ic) are given in 
the last row. 
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FIG. 1 : SAR is shown as a function of particle radius. Resonance 
can be observed at 9.21 nm with SAR = 6365 W/g. 

ation process exclusively did not explain the experimental 
data. Contrary, taking Neel relaxation alone, resulted in slight 
underestimation of measured values, the error was within 5%. 
It might worth be note, the absolute value of variation of the 
particle size from the mean value is nearly constant, which 
could be explained by having samples fabricated with the 
same technique. Contrary the samples of Mehdaoui et al. 
were prepared differently, and according to it, their standard 
deviations of particle sizes differ markedly. The asymmetry 
of the variation of particle size from the mean value can be 
explained by the asymmetry of the log-normal distribution. 

Wang et al. also estimated the optimum size for reaching 
the maximum of SAR being at r = 9.15nm, based on the 
linear response theory and Rosensweig's susceptibility. From 
our model, we got r = 9.21nm for the optimum size to get 
resonance, which size is in excellent agreement with the esti- 
mations given by Wang et al.. We note, that at the resonance 
the SAR is 6365W/g, see Figfl] The same curve can be seen 
on logarithmic scale, Figj2] to highlight the tendency for small 
nanoparticles. 

IV. SUMMARY 

The aim of this paper was to show the effect of size polydis- 
persity on the specific absorption loss in a magnetic nanopar- 
ticles containing media. For the first time, analytic calcula- 
tions were carried out instead of numerical simulations, using 
a Bloch-Bloembergen alike model, developed and described 
detailed in paper 02OI1 . As a reference, the paper of Mehdaoui 



et al. was used for direct comparison, including nanoparti- 
cles r > 7nm, the range of Stoner-Wohlfarth model beyond 
linear response theory. We found excellent agreement, see Ta- 
bled] Based on this agreement, we estimated the possible mean 
value and standard deviation of particle sizes for other experi- 
ments of Wang et al., Hergt et al. and Suto et al, see TablefTTTJ 
For the sake of comparison, the SAR values were calculated 
also without polydispersity, TableTTA and we found drastic im- 
provement when polydispersity was included. Our analytical 
approach predicted the resonance in the experiment of Wang 
et al.. After analyzing standard deviations from mean value 
of particle sizes, samples fabricated by different technologies 
might be distinguished, see experiments of Mehdaoui et al. 
and Wang et al. It can be concluded, the polydispersity is in- 
evitable in experiments and it should be taken into account in 
models for calculating SAR as well. The case of pure Brow- 
nian relaxation is excluded. Neel relaxation was found to be 
essential in relaxation processes for paramagnetic nanoparti- 
cles at hyperthermia. 
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FIG. 2: SAR is shown as a function of particle radius at logarithmic 
scale. Low SAR data are enhanced by logarithmic scale. Resonance 
can be observed at 9.21 nm with SAR = 6365 W/g. 
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